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Extracting B K* Form Factors from Data 

Christian Hambroclfl and Gudrun 

Institut fiir Physik, Technische Universitdt Dortmund, D-44221 Dortmund, Germany 

We extract ratios oi B ^ K* form factors at low hadronic recoil from recent data onB ^ K"" fi^ ii~ 
decays in a model-independent way. The presented method will improve in the future with further 
(angular) studies in semileptonic rare B-decays and advance our understanding of form factors, 
which are important inputs in precision tests of the Standard Model. 



Introduction. The experimental studies of flavor 
changing neutral current (FCNCs) decays of 6-flavored 
hadrons are advancing rapidly with the successful start 
of the Large Hadron Colliders' (LHC) 6-physics pro- 
gram. This generates further demand for good control 
of hadronic uncertainties to enable precision tests of the 
Standard Model (SM). 

A leading source of theoretical uncertainty in semilep- 
tonic FCNC decays B -> MU, i = e, /j.,r,u of a B me- 
son into a light meson M are the transition form factors. 
When the energy of the emitted meson is large in the 
B meson rest frame, the method of light cone QCD sum 
rules (LCSR) applies and allows to calculate the form fac- 
tors in this kinematical regime [1 . The region in which 
the emitted meson is softly recoiling against the decay- 
ing B meson is accessible to lattice QCD. Preliminary 
unquenched results are available [2 . As we rely on such 
predictions to probe elect roweak and flavor physics, it is 
important to obtain further independent information. 

The low recoil region allows for observables in which 
the form factors can be probed model-independently [T. 
Two of these short-distance free observables, the frac- 
tion of longitudinally polarized vector mesons Fl and 

(2) 

the transverse asymmetry Aj. |4j, have recently been 
measured by CDF [5 , BaBar [6] {Fl only) and LHCb f?! 
in the decay B K* ji^ ii~ . Here, we extract ratios of 
B iC* form factors from these data. Since the ratios 
are obtained at low recoil they can be directly compared 
to lattice predictions. It is the aim of this letter to demon- 
strate the performance of this new model-independent 
method from the first time available precision data. 

B K*i~^i~ at low recoil. The transversity am- 
plitudes ^^'Ho in 5 ^ K*l^l~ decays at low recoil 
factorize at lowest order [8j and share universal short- 
distance coefficients C^'^ 3 



i=±,||,0, 



(1) 



while the form factors fi are independent of the Wilson 
coefficients of the electroweak theory. A fourth transver- 
sity amplitude arises for finite lepton mass m^. However, 
at low recoil it is suppressed by m|/m^, where rriB de- 
notes mass of the 5-meson, and can be safely neglected. 
(We suppress the dependence on the dilepton invariant 
mass squared q'^ for brevity in this work.) 



It follows that observables of type {AfAf* ± 
)/(AMf ± ), where ij,k,l =±,||,0, 

probe long-distance physics only. This includes Fl and 
A'^\ given as 



E 
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X=L,Ri\Ao\'^ 



lL|2 



T |^L|2 + |^^|2 + |^L|2 + |^ii|2- 

Using Eq. ([l]) one obtains at fixed [3 , 

fi 



Fl = 
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The relevant B —5- K* g^-dependent form factors V,Ai 
and A2 are defined as usual 19 : 



2V 



{K*{k,e)\s^,b\B{p)) = 

rriB + rriK* 

{K*{k,e)\sj^j,b\B{p)) = 



re*Pp''k\ (7) 



2Ao m.K' + Ai {nfiB + niK* ) {g^ip - ) 
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where e denotes the polarization vector of the K* . 

It is clear from Eqs. Q and ([s]), that one can ob- 
tain ratios of form factors V/Ai and ^1/^42 by resort- 
ing to the observables Fl and Aj.. We comment on 
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the Umitations of this determination: The universality of 
Eq. ([T]) is broken by power corrections at order l/m^ to 
the improved Isgur-Wise form factor relations between 
the form factors of the dipole currents Ti 2,3 [9 and the 
form factors of the axial and axial- vector currents. How- 
ever, these corrections enter the transversity amplitudes 
with additional parametric suppression by the Wilson co- 



efficients Cn of the \AB\ 



\AS\ 



1 effective Hamil- 



tonian C^/Cg ~ (9(0.1) O [8]. Hence, we expect cor- 
rections at most at the level of a few percent, as for 
the l/m^-corrections to the operator product expansion 
(OPE) [To]. Further corrections to Eq. ([T]) arise in exten- 
sions of the SM with right-handed currents, which induce 

a finite at low q^. On the other hand, if SM-like 
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chiralities are present only, Aj. vanishes up to tiny cor- 
rections. The latter feature is in agreement with present 
data Oil], although within sizeable uncertainties. Im- 
proved data and global analyses [11, 12 including checks 
of the OPE [3 will exhibit whether such effects, which 
we neglect here, are contributing to AB = 1 transitions. 

Analytic form factors. For an analytic ansatz for 
the g^-shape of the form factors we employ the Series 
Expansion (SE) pISHlQ] . It is based on an expansion in 
z{t), defined as 



z{t) = z{t,to) = 



(9) 



and t being the analytic continuation of q'^ to the complex 
plane, with \z\ < 1. Here, t± = (m^ ± m^*)^ and < 
to < t+ for which we use to = t-^{l — ^/l — t-/tj^) [19l 
[^. 

We define the reduced transversity amplitudes fi = 
fi/M^ i =±, 11,0, where the g'^-dependent normalization 
factor M is given in Eq. (|6|. The SE can then be written 
as: 



m 



^a,,fcz^(t),(10) 



k 



with / = 1, 0, and m = 1, 1, for z =±, ||, 0, respectively. 
The Blaschke factor B{t) = z{t^m'^) parametrizes an 
off-shell pole associated with a meson with mass ttir; 
there can be several such factors if multiple resonances 
are present. The function ^/(t) is given by 



(t - 1+) 



(3+n)/2 



487rx/(n) {t^-t^y/^ 
to-t 



(11) 



with isospin factor = 2, and Wilson coefficients X/(^) 
with n being the number of subtractions in the renor- 
malization procedure (here n = 2), discussed in detail 
in [19l[20]. We take x/(2) = 1.2/(100mg). Note, that 
the ratios of the form factors, which are accessible with 
our analysis, are not sensitive to this input. 
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TABLE I: Recent high-^^ data from BaBar 6 , CDF [5 and 
LHCb with the statistical and systematic uncertainties 
added in quadrature. The maximum value in units of GeV^ 
equals X = 19 for LHCb and is the endpoint otherwise.^ Using 
A^^^ =2Ss/{1-Fl). 

We are constricted to 4 d.o.f., i.e., two bins in the high- 
q'^ region for each observable, see Table |lj Fortunately, we 
obtain a good fit already for the lowest order SE, where 



(t) = A(t, ml- ) ^y-z{t,0)^yz{t,t_) , 
/|| = a\\ K{t,ml^)^-z{t,{)) , 
/o(t) = ao A(t,mi+), 



with 



A(t,m^) 



M 



(12) 



(13) 



z(t,m2,)(/>^-^(t)' 

We take mi- = 5.42 GeV for the vector (±) and mi+ = 
5.83 GeV for the axial vector (||,0) transitions [21]. 

Fit and results. Even though Eqs. Q is exact for 
each (7^, the data on the other hand is binned. In the 
fit this is taken into account by replacing ff{q^) with 
j^.^^dq^ {2 pi{q^)ff{q^))^ where the universal function 
pi'= (|C^p + |C^p)/2 drops out of Eqs. Q in the limit 
of vanishing bin-size. To the accuracy discussed after 
Eq. (Isl) we have 



where 



■n^-^Cf{q^) 



\Ci,\\ (14) 



1 - 2^ In (/i/m^) - 1, and Cf 
denote the effective coefficients, see [3^ for details. Fig. [2] 
(left-hand plot) shows the slope of pi in the SM including 
next-to- leading order QCD corrections as used in the fit 
(solid curve) and with physics beyond the SM (shaded 
areas) . The latter is estimated by a variation of the Wil- 
son coefficients C7 from -0.4 to -0.3 and |Cio| from 2 to 
5, see [12]. Since the bin- averaged shift in the slope is at 
the percent- level, the uncertainties related to potential 
contributions from physics beyond the SM are negligible 
given the current accuracy of the data. 



The SE in Eq. (12) is truncated after the first order 
term, introducing three parameters ay and a^. Since 
the data constrains only ratios of the , we are sensitive 
to the modulus of the ratios of the ai only. We choose to 
fit and a^/a^ and obtain Xmin — ^-^ data, 

see Table [ij The error correlation of the fit is shown in 
Fig. [2] (right-hand plot), the fit to the data can be seen 
in Fig. [ij The best-fit results with Icr errors are 



<^||/<^± = 0.43 



+0.11 

-0.08' 



Oi{)/a±_ = 0.15 



+0.03 
-0.02 



(SE). (15) 
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FIG. 1: Experimental data by BaBar fBl (orange triangles), 
CDF [5 (blue circles) and LHCb 7 (black squares), see Table 
[l] versus the outcome of the fit at 68% CL (dark-green bands 
unbinned, dashed magenta box binned) 95% CL (light-green 
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bands unbinned) for Fl (top) and Aj. (bottom) 
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FIG. 2 
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Left: The function pi{q^) normalized to its value 
at ^min = 14 and 16 GeV^ in the SM (solid curve) and with 
non-SM effects (shaded region). Right: Error correlation of 
a\\/a± vs. ao/a±. The shown regions are 68% CL (dark- 
green), 95% CL (light-green) with the best-fit values given by 
the black square. 



The resulting form factor ratios V/Ai and ^1/^42 are pre- 
sented in Fig. [3j We further show the form factors them- 
selves, fixing a± = —0.08, i.e., V ^ to allow for a com- 
parison of the g'^-shapes. As can be seen, the sensitivity 
to A2 is very low towards the endpoint, caused by the A 
factor in the transversity amplitude /o, see Eq. ([5|. Also 
shown are LCSR predictions for < 14 GeV , where 
we assigned flat uncertainties as obtained at q'^ = 0, i.e., 
11%, 12% and 14% for V,Ai and A2, respectively [9]. 
Quenched [22 (set 3) and unquenched lattice results [2^ 
are shown as well after using the improved Isgur- Wise- 
relations to lowest order in 1/ rrib O |8] 



T,{q')=KV{q'), Uq') = ^A,{q'). 



(16) 



y(o)/Ai(o) 2.0 ±0.4 3.iig:? 1.4 ±0.2 1.5 ±0.9 
Ai(o)/A2(o) i.2±o.i i.iio.i 1.1 ±0.2 1.0 ±0.7 



TABLE II: Form factor ratios at q'^ 
and added in quadrature. 



0. * Errors symmetrized 



We add theory uncertainties in quadrature, thereby ig- 
noring possible cancellations of systematics in the ratios. 
Both lattice and LCSR predictions agree within 95%CL 
with the ratios obtained from the SE fit. 

Using the parametrization of the Simplified Series Ex- 
pansion (SSE) [18j, corresponding to the changes B{t) 



l-t/mji, <^/(£) ^ 1, ^/-z{t,0) Vi/mB, \/z{t,t-) 
in Eq. il2h, a similarly good fit ist obtained with 
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X: 



3.4 and the best-fit results 



= 0.44t^;J^, ao/a^ = 0.321^4 (SSE). (17) 



The form factors from the SSE are in good agreement 
with the SE ones and not shown. Noticeable differences 
arise only in the extrapolation to small g^, where the SSE 
results for V/Ai can be about 50% larger than the ones 
from SE, however the predictions still overlap within la. 

Form factor ratios at = along with LCSR predic- 
tions O |23] are compiled in Table [llj Both fits return 
values of V/Ai larger than the LCSR ones and larger 
than the symmetry-based prediction in the large energy 
limit V/Ai {rriB + rnK-Y/{m\ + ^1:* " El]^ 
z.e.,"l/(0)/Ai(0) = 1.3 up to l/m^ corrections [25^. A 
finer binning is expected to improve the control of the 
g'^-shape and the performance of the fit results at low q^ . 
We stress, however, that already with present data a con- 
sistent picture over the whole physical range is obtained. 

Conclusions. Determinations oi B ^ K* form fac- 
tor ratios from present data and theoretical predictions 
are mutually consistent, see Fig. |3] The currently still 

— (2) 

early data may shift towards a larger Aj.^ at low recoil 
and lead to a correspondingly lower value of V/Ai, im- 
proving the agreement with the preliminary unquenched 
lattice determination. However, once data have become 
sufficiently precise allowing for consistency checks [iQl - 
12 , the low recoil predictions constitute a powerful test 
for the lattice and its advances [26 . 

The method of form factor extractions from 5-decays 
into a fully reconstructable final state will greatly benefit 
in the near-term future from the high-luminosity flavor 
searches at the LHC; it is applicable to other decays in- 
cluding Bs ^ (j) Both better statistics and the 
assessment of further short-distance free observables [3] 
from a full angular analysis [27 will improve the control 
of the low recoil region and of the form factors. 
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FIG. 3: The form factor ratios V/Ai and A1/A2 from fit to data 0{7] (upper row). In the lower row the form factors are 
shown using V as input with a± — —0.08. The best fit result with 68% CL (95% CL) bands are given by the solid black 
line and the shaded dark-green (light-green) region, respectively. The solid blue line and the hatched region denote LCSR 
results 9 . Lattice results from ^ and [22], using the Isgur-Wise relations Eq. (16), are shown as red circles and orange boxes, 
respectively. 
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